Quasiperiodic fields and Bose-Einstein condensation 
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Abstract 



We construct a partition function for fields obeying a quasiperiodic boundary 
condition at finite temperature, i/}(0;x) = e ie ^(P;x), which interpolate continously 
that ones corresponding to bosons and fermions and discuss the possibility of con- 
Q densation for these fields. 

<N 

oo 



PACS: 11.10. Wx; 05.30.-d; 03.65.Db. 

Keywords: Bose-Einstein condensation; partition functions. 



*e-mail: pborges@if.ufrj.br; 
'''e-mail: boschi@ctp.mit.edu; 
' : e-mail: farina@if.ufrj.br. 



The concept of condensation was proposed by Einstein , after the fundamental work 
of Bose [0 on the statistics of integer spin particles, and its first realization was given 
by London in a model for the Helium superfluidity [?]. More recently, the interest for 
this phenomenon has increased because of the experimental results on atomic traps and 
the direct observation of Bose-Einstein condensation (BEC) on gases |4| signalizing for a 
large number of applications. 

Long ago, it was shown for massive bosons in the absence of external fields that BEC 
can only occur if N > 3, where N is the space dimension. For massless bosons this 
condition reduces to iV > 2, although in this case the critical temperature is infinite. 
These results were first shown in a nonrelativistic context and later generalized the 
relativistic case ||||. Very recently, the extension of these results to curved space-times 
and the inclusion of non-uniform magnetic fields in the relativistic context were also done 
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Here, we present a relativistic quantum field theory approach to BEC from which we 
obtain nonrelativistic results as a particular limit. We start from a partition function 
identified with a determinant for fields which satisfy a quasiperiodic boundary condition, 
i.e., a condition which interpolates continuously the periodic (bosonic) and antiperiodic 
(fermionic) cases at finite temperature. This kind of field is usually considered in lower 
dimensional systems as in the anyon theory in 2+1 dimensions [O, |I2|. Recently, we have 



shown |T3[ that the present construction leads to a partition function which is equivalent 



to the usual ones for anyons, at least for the second virial coefficient [T%\. Further within 
our approach, we also obtain a general condensation condition for fields which interpolate 
continuosly between bosons and fermions and discuss the role of the dimension where the 
condensation occurs. 

Let us start by writting the partition function for a general system described by a 
Hamiltonian H with chemical potential fi as 

Z = exp{-/3fi} = tr e ~^ H+ ^\ (1) 

where Q is the free energy and M is the conserved charge (conserved quantum number) 
of the system. As it is well known, for relativistic charged massive bosonic fields one can 
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express this partition function as a determinant || |15[, namely, 

-l 



Z 



Bosons 



det(-D 2 + M 



(2) 



where D 2 is the square of the covariant derivative, including the chemical potential, D v = 
(do + ifj,, &i) and M is the mass of the field. This prescription for the inclusion of the 
chemical potential /i coincides with the one given by eq. ([TJ) except for the introduction 
of an imaginary part to 2, which is not physically relevant. So unless otherwise specified, 
we are considering that the real part of Z has been taken. Note also that, the chemical 
potential appears squared in the Klein-Gordon (KG) operator, —D 2 + M 2 . This is a 
consequence of representing the partition function in phase space and then integrating 
over momenta (see || 16] for details). The label P means that the eigenvalues of this 



operator are subjected to periodic boundary conditions and hence they are given by 

Xnk = K + ifi) 2 + k 2 + M 2 , (3) 

where u n = 2mr//3, with n G 2, are the Matsubara frequencies [|I7| for bosonic fields and 
k G IR\ 

Then, we recall the partition function for a free realtivistic fermionic gas. As for the 
case of bosons, one can write this partition function as a determinant || |l^] 



Z 



Fermions 



det D (ip-M)\ A 
det(-D 2 + M 2 ) 



(4) 



where deto means the calculation over Dirac indices and the subscript A means that the 
eigenvalues of each operator are computed with antiperiodic boundary conditions. Hence, 
for the operator — D 2 + M 2 , the eigenvalues are given by (|3]), but now the Matsubara 
frequencies are u) n = (2n//3)(n + 1/2) with n£2. 

Now, with the aid of a convenient choice of parameters we are going to use an 
approach which describe both the bosonic and fermionic partition functions in an unified 
way. Consider then the generalized determinant 



Z 



(5) 



det(-£> 2 + M 2 

where the subscript 9 means that the quasiperiodic finite temperature condition, i.e., 

^(0;f) = e*V(/3;x), (6) 
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is assumed and we introduced a new parameter a to be able to reproduce correctly the 
bosonic and fermionic particular cases. Observe that when 9 = and a = — 1 we obtain 
the bosonic partition function, while for 9 = it and a = +1 we have the fermionic one. 
The eigenvalues in this case are 

2 



A 



2nix 9 
+ ifi + - 



+ k 2 + M 2 . (7) 



P ^ P. 

The determinant (|5|) is a generalization of its quantum mechanical (0 + 1) dimensional 
version (with \x = 0) which has been calculated using Green functions [I8| and the zeta 
function method fll9|| . 

Notice that the interpolating parameter 9 plays the role of the zeroth component of 
a topological (temperature-dependent) gauge field = (A (P),0), since the ^-boundary 
condition implied that u n — > u n + 9/(3, which can be viewed as a shift in the time 
derivative operator, do — > do + i9/ [3. The discussion of the bosonic (0) and fermionic (|J) 
determinants with a constant topological gauge field, but without relating it to boundary 
conditions or interpolating statistics, has been presented before ||16|| . It is important to 
mention that the parameter 9 or equivalently the Ao{(3) field can not be gauged away 
because of the nontrivial boundary conditions introduced by the time compactification 
to the interval [0,/?]. An interesting point that has not been observed before is that the 
coupling between bosons with the topological field A (j3) suppresses the Bose- Einstein 
condensation, unless for very specific values of 9, as we are going to show. 

It can be shown that the free energy corresponding to the partition function ([5]), is 
given by p|: 

n((3,n) = -iln-Z = -^lndet^. 

= AaV [^) E cos M) cosh (^)(-) K h(N+1) (n(3M) , (8) 

where K u (x) is the modified Bessel function of the second kind. Note that in the above 
formula, only the real part of the free energy f2(/3,/i) is expressed, according to the 
prescription of introducing the chemical potential as an imaginary time-component gauge 
potential H. If we particularize the parameters a and 9 to the bosonic (a = — 1 and 
9 = 0) and fermionic (cr = +1 and 9 = tt) cases we shall find precisely the results known 
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in the literature [0, |2(| [16|]. In a recent work |I3 |, we have shown that this partition 
function with arbitrary 9 leads, in the nonrelativistic limit, to a second virial coefficient 



which is equivalent to the well known result of Arovas et. al JT4J. We have also calculated 
higher order virial coefficients showing good agreement with anyon perturbation theory 
and numerical calculations [pl |. There, we have shown that, for statistical consistency we 
have to impose that the parameters 9 and a are related by 



a 



a{9) 



cos" 1 9 



(9) 



To see explicitly when the condensation occurs, let us rewrite the free energy (||) as |L6| 



_TVa(9) 

TVa(9) 
2(2tt) n 



d N k K {in 



1 - e w e 



+ -»■ -»)} 



d N k{\n[l + e~ 2(3{u) -' M) -2cos9e-^-ri] + -»■ -//)} (10) 



where u = y k 2 + M 2 . Further, the charge density is given by: 

1 d , _ 1 



P = 



InZ 



IT 



j3V d[i " V \d[i J VT 

Remind that it goes to infinity when particles condense. From this charge density we can 
write the general ^-dependent distribution function 

1 



fe(J3,±fi) 



1 

~2cos# 



cos 9 - e~^ T ^ 
cosh (5{uj =f A*) — cos # 



(12) 



The condensation condition comes from the divergence in the free energy or equivalently in 
the corresponding charge density which is also evident in the above distribution function. 
As the condensation is related to the zero momenta (k = 0) state, this implies that u — ► M 
and in the general case, with arbitrary 9 we have that this condition is given by 



or simply 



1 + e-W-ri - 2 cos 9 e-PW^ = 



cosh/3(/x — M) = cos 9 



(13) 



(14) 
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which can only be satisfied for finite temperatures if 9 = 2mr (with n integer) and fi = M, 
simultaneously. This means that the condensation condition in the general case is exactly 
the same as in the usual BEC. 

Then, we conclude that equation (|10|), or equivalently (|j), is well defined^ for every 
value of /i, except when 9 = 2mx for which [i = M gives the well known condensation 
condition. So, we can infer that no condensation occurs for any 9 ^ 2nir (n integer). 
However, this does not forbid the condensation of fermions or ^-fields. 

Consider, for instance, that the free #-field obeying the Klein-Gordon equation is now 
interacting with another topological field A' Q , different from the A Q introduced above to 
transmutate the statistics of the original bosonic or fermionic field. In this case, the 
parameter defining the statistics and condensation condition for the ^-field together with 
the topological interaction can be written as B = 9 + 9', where 9' is associated with the A' 
field . Then, the statistical properties in this situation are determined by the parameter 
O. This means that apart from the usual BEC, fermions coupled to topological fields 
with 9' = nn so that the total phase is G = 2nn may possibly condense. Analogously, 
any #-field coupled to a topological field with 9' = 9 — 2mr, such that the total phase is 
G = 2nir may also condense. 

To obtain the critical temperatures, densities and dimensions for the condensate from 
the above discussion we go back to the free energy given by Eq. (||). Note that this 
equation contains all the relevant information about our system. However, we are mainly 
interested in the ultra and nonrelativistic limiting cases. 

Let us first consider the ultrarelativistic case. In this situation where the energies 
involved are much higher than the mass scale (the mass M of the particles) we should 
take the limit (3M « 1, which is easily recognized as a high temperature limit. In this 
regime we can use the following property of the Bessel functions 

\xmKJx) = Y(v)2 u ^x~ u ; (u > 0) (15) 

x— >0 

1 Note that from ([h]) or ( |l2| ) one can see that there is another possibility of infinite charge density 
which comes from the choice 6 = (2n + 1)tt/2, since a — > ±oo, because of relation (0). However, this 
divergence does not depend on the other physical parameters of the system and then this sigularity 
opposed to the BEC case, seems to be unphysical. 
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so that the free energy (f|) reduces to (9 = 2nn, a{6) = — 1) 

2V N 4- 1 + °° 1 

= - ( ^ r+1 r(^-) £ cosh(n/^)(-)" +1 . (ie) 

Using now Eq. (|llj) we find that 

2 / l\ N N 4- 1 + °° 1 
^^)=(^pTT^j r(^-)Esinh(n/3 / i)(-) Jv . (17) 

As we are dealing with the high temperature limit, (3M « 1, the condensation condition 
fi —> M also implies (3[i « 1, so that 

p=Jr^n^)aN-i)T«-\ (is) 

where ((s) = J2^=i n ~ s is the usual Riemann zeta function. Using the fact that the critical 
temperature of the condensate is reached when /i = M, in N = 3 space dimensions, we 
have: 



which coincides with the ulatrarelativistic BEC results known in the literature [0-|lC]] • 

Now, let us study the situation for the nonrelativistc limit which here means to take 
the low temperature limit, 13 M » 1. Using the asymptotic expansion for the Bessel 
function 

K v {x) ~ yje-*, (20) 

valid when |x| >> 1 and —tt/2 < argx < tt/2, and taking f3M » 1 in the free energy(|8|), 
we have (6 = 2nix, a{9) = —1): 

(„ A' 
f-) ^sh(n^)(-Y^e-^ M . (21) 

In this regime the charge density reads 



N 



(J\/f \ 2 +00 1 
—J Esmh(n/3/i)(-)fe-^ M . (22) 

Exactly as in the ultrarelativistic case, here the critical point is reached when /i — > M. 
So, in the nonrelativistic limit, the condesation condition implies that (3fi » 1 and then 
the charge density reduces to 

JV 

^,)=g) T f(i)^-»). (23) 



At the critical point fi = M we get simply 




(24) 



which for N = 3 space dimensions leads to the critical temperature 



27T (_p_ 



) 



2/3 



T, 



c 



M {((3/2) 



(25) 



in agreement with the nonrelativistic BEC results known in the literature PHTOf. 



ultrarelativistic and nonrelativistic cases one can see that the condensate is not defined 
in two space dimensions for both cases, since the Riemann zeta function ((s) has a pole 
at s — 1. Hence, we find that the condensation for free massive fields (both in the ultra 
or nonrelativistic limits) occurs for space dimensions N > 2, which is in accordance with 
the BEC results for the massive case in the absence of external fields |)|-|l(J . 

We presented here a relativistic partition function from which we described both the 
nonrelativistic and ultrarelativistic condensation. With our partition function we were 
also able to describe fermions and bosons at the same time by varying the parameters 9 
and a. From this partition function we found the general condensation condition of the 
fields, i.e., they should strictly obey periodic boundary conditions at finite temperature. 
This condition can be achieved by any field including fermions but conditioned to the fact 
that they must be coupled to topological fields which result in a total phase G = 2mr. 
On the other side, bosons coupled to those fields in such a way that ^ 2nn do not 
condense. 

We also described here a kind of finite temperature bosonization (or fermionization) by 
including the arbitrary 9 in the partition function. We have shown that no condensation 
occurs for fields with 9 ^ 2nii. However, a fermionic field interacting with a topological 
field A = tt//3 will also condense, so this should correspond to a bosonization process. 
Reversely, adding such a field to a boson we fermionize it once they will exclude rather 
than condense. Note that this picture is valid only at finite temperatures since at zero 
temperature the topological field A = 9/(3 vanishes. 

Our results also suggest that pairs of fermions can also condense as they would satisfy 
the condition O = 2nir, as well as arbitrary combinations of ^-fields as long as that 



Looking at the Eqs. (|IS|) and (^) for the charge density in N space dimensions for the 
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condition is preserved. However, to have a clear picture of such situations one needs to 
consider some kind of interaction between these fields which was not considered in our 
approach. This is presently under investigation and will be reported elsewhere. 
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